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Abstract 

This paper provides the quantum treatment of the relational quadrilateral. The underlying reduced configuration 
spaces are CP 2 and the cone over this, C(CP 2 ). We consider exact free and isotropic HO potential cases and pertur- 
bations about these. Moreover, our purely relational kinematical quantization is distinct from the usual one for CP 2 , 
which turns out to carry absolutist connotations instead. Thus this paper is the first to note absolute-versus-relational 
motion distinctions at the kinematical rather than dynamical level, and also an example of value to the discussion 
of kinematical quantization along the lines of Isham 1984. This treatment of the relational quadrilateral is the first 
relational QM with very new mathematics for a finite QM model, and is far more typical of the general quantum 
relational iV-a-gon than the previously-studied case of the relational triangle. We consider useful integrals as regards 
perturbation theory and the peaking interpretation of quantum cosmology. We subsequently consider problem of time 
applications of this: quantum Kuchaf beables, the Machian version of the semiclassical approach and the timeless 
naive Schrodinger interpretation. These go toward extending the combined Machian semiclassical-Histories-Timeless 
Approach of [1] to the case of the quadrilateral, which will be treated in subsequent papers. 
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1 Introduction 



Following on from Paper I's [1] classical work on the quadrilateralland relational particle model (RPM) [2, 3], 1 the 
present paper considers quantum relational quadrilateralland. Quantum RPM's were first considered by Barbour, Smolin 
and Rovelli [4, 5], though practical progress with solving concrete examples of these was hindered until [6, 7, 8]; see 
[9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 1, 21, 22, 23, 24, 25, 26, 27] for subsequent development. This sudden progress 
was triggered by Barbour formulating the pure-shape RPM [3, 28] (which proved to be easier to solve) and by the sequence 
of Keys and profitable interdisciplinary observations given in Paper I (for which [29, 30, 31, 32, 33, 34, 35, 36, 37] are 
particular antecedents from Molecular Physics, Celestial Mechanics, Geometrical Methods and Shape Statistics). 

RPM's have been valued as toy models [38, 39, 40, 19, 41, 42, 22, 23, 25, 24, 26] for the Problem of Time (PoT) 
[38, 43, 23] in Quantum Gravity via the analogies exposited in [44, 25]. They have been used e.g. for study of Timeless 
Approaches [45, 39, 11, 14, 15, 17, 18, 46, 25], the Semiclassical Approach [9, 17, 18, 20, 25, 27, 47] in close parallel to 
Halliwell-Hawking's approach [48] to Quantum Cosmology, and Histories Theory and its combination with the previous 
two [21, 25]. 

In more detail, as argued in [23], the Temporal Relationalism and Configurational Relationalism which RPM's are 
constructed to embody are 2 of the 8 PoT facets. Which classical PoT aspects are present for RPM's and how to resolve 
them was covered in Paper I. At the quantum level, if Configurational Relationalism was resolved at the classical level, 
it stays resolved. Temporal Relationalism now resurfaces as the Frozen Formalism Problem; various strategies for this 
were laid out in Paper I; the main one followed in the present paper is the Semiclassical Approach [25, 27]. The Problem 
of Beables was resolved at the level of classical Kuchaf beables in Paper I. The present paper considers the quantum 
counterpart of this. As regards other PoT facets, the following statement holds at both the classical and quantum levels. 
RPM's are conceptually free of the Foliation-Dependence and Spacetime Reconstruction Problems, and can be cast as 
free from the Constraint Closure Problem too. That covers the 6/8ths of the PoT required to have a local resolution, 
with two caveats that require further papers (III and IV). 

1) The Semiclassical Approach requires support from other approaches, in particular Histories Theory and Timeless 
Records [49, 50]. 

2) Quantum Dirac beables have not yet been considered for quadrilateralland; the current program's [25] way of addressing 
these uses the machinery of 1) and thus must await treatment of that. 

RPM's and similar are also often used to motivate the Linking Theory Approach to Shape Dynamics [41, 51, 26], whilst 
other papers found uses in investigating such as whole-universe path integral approaches [42] , geometrical quantization [5] , 
operator-ordering in Quantum Cosmology [52, 25], and an investigation of an alternative anomaly-based emergent-time 
mechanism [24]. 

Quantum RPM's hitherto studied are scaled 3-stop metroland [7, 9, 16, 20], pure-shape 4-stop metroland [14] pure- 
shape triangleland QM [13, 15], scaled 4- and TV-stop metroland [17], and scaled triangleland [18, 20, 21] . Sec Part III of 
[25] for a review of these. Quadrilateralland QM (pure-shape and scaled) is then logically the next step for this program 
and that taken here in the present article. Compared to the preceding list, there is now nontrivial CP fe ^ SU(k + 1) 
mathematics to contend with. The A-a-gon is unlikely to prove much harder than the quadrilateral, whereas the triangle 
is exceptionally simpler by benefitting in non-generalizable ways because its CP 1 = S 2 allows for extra techniques. 

Some interdisciplinary comments on the present Paper are as follows. [54, 53] consider the atom in N-d (in the sense of a 
1/r potential in dimension N) . [53] considered the Stark effect not only for A-dimcnsional atoms but for TV-dimensional 
rotors as well. Both of these are maximally symmetric problems (on R p and S p , each of which possess p{p + l}/2 
Killing vectors). MacFarlane's work [37] and the current paper can then be viewed as an extension of this work for the 
next most symmetric case of CP 2 that exists for shape space dimension q = 4 ( this has 8 Killing vectors rather than the 
maximal 10). Our paper's useful integrals for QM on CP 2 further extend MacFarlane's work to perturbations about the 
free case. See [55] for other literature concerning HO's on CP W , though we do not know of any previous literature that 
covers the CP 2 counterpart of the Stark effect. It is nontrivial as a robustness test of the atom, in that it unveils a number 
of fortunate occurrences for the standard orbitals and Stark effect that end upon passing from maximal to the next most 
maximal symmetry. See the Conclusion for examples of this. One part of the interpretation of QM of a quadrilateralland 
involves an application of Paper I's complex-pro jective chopping board counterpart of Kendall's spherical blackboard from 
Shape Statistics [33] . We shall also see that this QM is a cross between the Periodic Table and Gell-Mann's eightfold way 
from Particle Physics, in a sense made precise in Sees 5 and 6. We shed light on how quadrilateralland's HO-type systems 
are far more like triangleland 's than 4-stop metroland's at the quantum level, despite their greater classical similarity 
with 4-stop metroland indicated in Sec 1.16. This has further relevance as regards 'triangleland within quadrilateralland' 
robustness tests paralleling Kuchaf and Ryan's work [56] in minisuperspace Quantum Cosmology. 

An outline of this Paper is as follows. In Sec 2, we consider kinematical quantization for pure-shape quadrilateralland. 
It is an interesting example as regards Isham 1984 kinematical quantization [57] and as regards how the absolute versus 
relational motion debate already shows up at the level of kinematical quantization. (One of us previously pointed out 
distinctions of this type at the subsequent level of the wave equations themselves [25].) In Sec 3, we construct the 

Refer to Paper I for details of motivation for RPM's in general and quadrilateralland in particular, for notation and for cqua- 
tion/section/figure references that begin with 'I.' 
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conformal-ordered TISE for pure-shape quadrilateralland, which we separate for the free case in Sec 4 in Gibbons-Pope 
type coordinates [32, 37]. Sec 5 covers the energies, quantum numbers and wavefunctions for this problem and Sec 6 
describes the ground state and first few excited states. Sec 7 considers scaled quadrilateralland, in particular for isotropic 
HO's. Sec 8 constructs useful integrals out of the scaled isotropic HO and pure-shape free wavefunctions (these generalize 
the integrals used in e.g. the study of the Stark effect). These are then used in Sec 9 for peak and spread analysis 
('Peaking Interpretation of Quantum Cosmology', though we show this has a well-known counterpart in Atomic Physics) 
and Sec 10 for time-dependent perturbation theory in the Semiclassical Quantum Cosmology analogue model context. 
Sec 11 finishes Paper Fs consideration of Naive Schrodinger Interpretation questions for quadrilateralland; this approach 
is a prequel to the Machian version of Halliwell's combined approach's use of regions in Paper IV. We conclude in Sec 12, 
including a sketch of extensions to the general iV-a-gonland. 



2 Kinematical quantization of quadrilateralland 

In general, one has to make a choice [57] of a preferred subalgebra of functions of one's configurations Q c and momenta Pc 
that are the ones to be promoted to QM operators. Some context for this is that the Groenewold-van Hove phenomenon 
[58] precludes simultaneous promotion of all classical quantities to quantum operators. There are also global considerations 
[57] by which a model's quantum commutator algebra is not in general isomorphic to that problem's classical Poisson 
bracket algebra. 

Key 18 The RPM program lies within Isham's [57] Q/G example for G a subgroup of Q. Then the relevant spaces involved 
in kinematical quantization can be decomposed as semisimple products V*(Q) © G can (Q). Here, G can (Q) is the canonical 
group and V* is the dual of a linear space V that is natural due to its carrying a linear representation of Q that realizes 
the Q orbits. Mackey Theory is then a powerful tool for finding the representations of such semidirect product algebras 
[57]. Furthermore, V* = V for finite examples and G can (Q) = Isom(Q) for all 1- and 2-d RPM's. 

Example 0) For absolute W, the canonical group is Isom(RP) = Eucl(p) = Tr(p) © Rot(p) = W © SO(p). 

Then an appropriate linear space is R p , so, overall, one has R n ® R n ® SO(n). These are the x z , their conjugates 
the pi and the corresponding angular momenta Li = e^x^pk, so this case is both physically and mathematically very 
familiar. 

Example 1) For scaled TV-stop metroland, 1Z(N, 1) = R™, so the outcome is mathematically the same as above. However, 
physically the roles of the objects involved are relative Jacobi separations p l , their conjugates 7Tj and relative dilational 
momenta 2?il r (Sec 1. 10) for T running over 50(n)'s 1 to n{n — l}/2 indices. 

Example 2) For scalefree TV-stop metroland, S(N, 1) = S n_1 , for which the canonical group is Isom(§" _1 ) = Rot(n) = 
SO{n). Then an appropriate linear space is K™. Now the objects in question arc the 2?ilr again, alongside the n l 
that square to 1 so as to provide the on-S™ -1 condition. These unit Cartesian vectors in configuration space are most 
conveniently expressed in ultrasphcrical coordinates (since the 2?ilr are). 

Example 3) For pure-shape iV-a-gonland's S(N, 2) = CP™ -1 shape space, the canonical group is G can (S(7V, 2)) = 
Isom(CP™ _1 ) = SU(n)/Z n Moreover, this shape space can also be written as SU(n)/U(n - 1); thus it is also a subcase 
of the general form in Isham's example above. 

Then one possible kinematical quantization involves V© G can = SU(n)/Z n © R 2n , for R 2n better thought of as C" [57]. 

Note however that triangleland admits a distinct kinematical quantization, V©G can = 50(3) ®M 3 with the R 3 made 
up of the Dragt coordinates ([30], 1.29) Dra r and the SO(3) of mixed relative angular momentum and relative dilational 
momentum quantities as per Sec 1. 10. Moreover, this alternative i) does not involve postulating objective existence to 
absolute entities (present among the C 2 of relative Jacobi vectors) and ii) is a more minimal realization (3-d to 4-d). 

Moreover, generalizing the latter relational kinematical quantization of the triangle to the quadrilateral is not partic- 
ularly obvious. Progress can be made via noting that R 3 is also IHP(C 2 , 2) - i.e. the space of irreducible homogeneous 
polynomials of degree 2 (Sec 1.8) - via the 3-vector to Pauli matrix map that rests on the well-known accidental relation 
between SU(2) and 50(3). IHP(C",2) then continues to be available for general-n CP™ -1 kinematical quantization. In 
the quadrilateralland case, this space is composed of the 8 independent shape quantities of Sec 1.8. All in all, we have 
the kinematical quantization V © G can = 5£/(3) ® IHP(C™, 2) This is not the minimal-sized space for any N > 3 since 
2n < n 2 — 1 for all integer n > 2. Nevertheless, minimality is a guideline and not an obligation, and argument i) continues 
to stand. 

Next, the quadrilateralland isometry generators are Tp = {y,Z3,Z±,M±,V±}> among which those with particularly neat 
expressions are 

~ . , d cos 6 [ d n d 1 ~ , d sm.6 f d n d 1 ,_. 
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Finally, 

fa ( 1 8 . n d 1 f d 2 n „ d d d 2 I I 

we give this in the operator-ordering that is relevant to this paper's time-independent Schrodinger equation to be in terms 
of the Laplacian (see the next Section) . These can then be paired with the Gibbons-Pope type coordinate expressions for 
the shape quantities s r of Sec 1.6 so as to demonstrate closure and evaluate the commutators. 

The kinematical quantization of GR-as-geometrodynamics itself involves C°°(£, M(3, K)) © C°°(S, GL + (3, K)) [57], where 
the latter factor is closely associated with the mathematical identity of Riem(S), GL stands for 'general linear' and M(3, R) 
are real 3x3 matrices. 



Key 19 One can furthermore consider the above choice of kinematical quantization as also a selection of relational beables: 

a subset {K&} of the classical Kuchaf beables K — F[x, P,4>, ipiPx'PP'P&'Pi'] then promoted to the quantum level {K/\\ 
such that 

0) They obey [K, Linz] = which is trivially the case here by prior classical reduction. 

1) They cover all the relational information. 

2) They obey suitable continuity conditions. 

3) They themselves algebraically close under the commutation relation. 

4) They are allowed some redundancy (meaning more relational functions than there are independent pieces of relational 
information - not to be confused with including unphysical/gauge/non-relational information). 

Note that 2) to 4) are e.g. already evident in using sin<^> and cos</> for the circle/3-stop metroland. 

Then a candidate for the algebra of Kuchaf beables for quadrilatcralland is that it is the same as the kinematical 
quantization algebra of the T T and s r . 



3 Time-independent Schrodinger equation (TISE) for quadrilateralland 

Operator-ordering is more of a problem in Quantum Cosmology than in Atomic Physics due to paucity of observations. 
Can theory alone determine operator ordering? 

DcWitt [59] considered elevating the classical coordinatization-independence of configuration space to additionally 
hold at the quantum level. This suggests that the classical kinetic term N ab (Q c )P^Pq be promoted to the quantum-level 
Laplacian operator ordering 

A:= - 1 J^-/,/m(Q c )N ab (Q c )^1 (4) 

(this is also advocated in e.g. [60]). Moreover, this is not a unique implementation of DeWitt's criterion since one can 
include a Ricci scalar curvature term so as to have, for any £ <E R, 2 the ^-operator ordering 

A« := A -£Ric(Q;M] . (5) 

Among these £-orderings, there is [61] a unique configuration space dimension q-dependent (q > 1) conformally-invariant 
choice of operator-ordering ([62, 63]), 

A c := A - ? Ric(Q; M] := A - ^— ^Ric(Q; M] . (6) 
This furthermore requires that itself transforms in general tensorially under conformal transformations [61], 

* — > $ = n {2 - q}/2 ^ . (7) 

The TISE following from the above family of orderings is then (Euni denotes the total energy of the model universe) 

= E Uni * => A c * = 2{V - E Uni }*/fi 2 . (8) 
What is the underlying conformal invariance in question? [E.g. it does not act on space itself.] 



Misner's identification [62] is that it is the underlying conformal covariance of the Hamiltonian constraint under scaling 
transformations. 

H = Q — > % := Q- 2 n = . (9) 

2 Even then, an underlying simplicity here is that the above is the extent of the ambiguity only if one excludes more complicated curvature 
scalars e.g. by stipulating no higher-order derivatives nor higher-degree polynomials in the derivatives. 
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This can be generalized to conformal covariance of other quadratic constraints such as the energy constraint £ or the r- 
formulation counterpart £ r . Key 20 E.A. 's identification [52], on the other hand, goes one level deeper to the consideration 
of actions. It then so happens that it is the conformal invariance 

ds 2 -> ds 2 = tfds 2 , E-V^{E Uni -V} = {E Uni -V}/0 2 (10) 

of relational product actions [c.f. (1.1) and (1.17) for examples]. I.e. it is a conformal invariance of the kinetic arc element 
ds alongside a compensatory conformal invariance in the potential factor W, reflecting that the combination actually 
present in the action, ds, is not physically meaningfully factorizable. One then recovers Misner's conformal covariance for 
the purely-quadratic constraint one's theory possesses as a primary constraint due to its relational product form. 

For the physical quantities to be invariant, the inner product in this convention is to have a weight function to scaling as 
(see e.g. [13]) 

UJ > LO = Q~ 2 Ld, (11) 

so that 

J DQ?[i 2 w = J BQQ g V* 1 y 2 ti 2{2 ~ q}/2 un~ 2 = J OQtfJ^w . (12) 

Example 1) For triangleland and 4-stop metroland, q = 2 so £ c = and conformal ordering = Laplacian ordering. 
Example 2) For the general CP" -1 , q = 2{n-l} and Ric = 4n{n-l}. Thus A c CP „-i = A CP „-i-2n{n-l}{n-2}/{2n-3}. 
In particular, then, for quadrilateralland A£, p2 = A CP 2 — 4. 

The corresponding TISE's are then as follows. 
For 4-stop metroland in terms of 2?Tot, 

P Tot * = A §2 * = 2{V - E}*/fi 2 . (13) 

For iV-a-gonland in general terms, 

A CP „-i4< - 2n{n - \}{n - 2}*/{2n - 3} = 2{V - E}^/h 2 . (14) 
Then specializing and further specifying for the triangle in terms of <Sr t) 

5r t* = A §2 * = {V - E}*/2fi 2 , (15) 
whilst doing so for the quadrilateral in terms now of two distinct 'felt charges' associated with I Tot and y gives 

1 A ( sin 3 XC0Sx ^\ _ ?™* ^* + Zn{n-l}{n-2} ^ = 2{E Vni - V} ^ 

sin 3 x cos x dx \ d\ J sin 2 x cos 2 x 2n - 3 h 2 

The GR-as-geometrodynamics counterpart of this is the Wheeler-DeWitt equation 

:= -fi 2 ' j-^^ jVMN^^} - ^MlAlVlj _ Vhftic(^; VI* + 2^A* = , (17) 

coupled to a QM momentum constraint in general. ' ' here indicates regularization, well-definedness and operator-ordering 
issues. The minisuperspace version has partial rather than functional derivatives, a dimension-dependent operator-ordering 
coefficient for RicM rather than the infinite-dimensional limit and no QM momentum constraint. 



4 Quadrilateralland QM separates in Gibbons— Pope type coordinates 

The hydrogen atom's TISE separates in both spherical and parabolic coordinates [64] and that of the isotropic HO in both 
Cartesian and spherical coordinates [66]. However, it is generally regarded as quite good fortune to be able to separate a 
quantum problem at all. How do RPM's fare? 

The free and isotropic-HO type problems, A^-stop metroland [14, 17] and triangleland [13, 15, 18] are separable in (ul- 
tra)spherical coordinates. iV-stop metroland is also separable for these and the diagonal anisotropic HO type problem 
in the Cartesian coordinates that physically represent the relative Jacobi inter-particle cluster separations [8, 25]. Scaled 
triangleland [18] is also separable for the above potentials and for the diagonal anisotropic HO in parabolic coordinates, 
that here physically signify a split into subsystems ('base and median'). For the quadrilateral, we have specifically consid- 
ered [37, 1] the Gibbons-Pope type coordinates as best-possible analogues of the (ultra) spherical coordinates, and indeed 
these do not disappoint when it comes to separability of the free (or isotropic-HO type) TISE: Key 21. Quadrilateralland 
separates in Gibbons-Pope type coordinates. These coordinates are SU(2) x [/(l)-adapted, and that part separates out 
as a package [using *(x,/3,0,^>) = A(ip,<j), (3)R(x)}- 
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The free problems for pure-shape 4-stop metroland and triangleland then form the even more standard package solved by 
the spherical harmonics (themselves separable into an SHM part and an associated Legendre part). Moreover, pinning 
physical interpretation on these, they come as, firstly, 

f> Tot *(0, 4>) = D{D + 4>) (pure-shape 4-stop metroland) (18) 

for D total relative dilational momentum quantum number, with second quantum number d the relative dilational mo- 
mentum between the base and median subsystems, featuring as the eigenvalue in T> 3 ^ = d<J>. Secondly, 

5 T ot*(e, $) = S{S + l}*(e, $) (pure-shape triangleland) (19) 

for S the total shape momentum quantum number: mixed relative angular and relative dilational momentum [15, 25]). 
Here, the second quantum number either is j the pure relative angular momentum component of the base relative to the 
median in the DES basis or an also mixed relative angular momentum and relative dilational momentum shape quantum 
number s in the EDS basis. This features in the eigenvalue problem = j \I/ or s^f (which depends on the meaning of the 
principal '3' axis in each coordinate basis). These are clearly analogues of the angular momentum eigenvalue equations 
for the rotor and atom. For quadrilateralland, 

l Tol A^,(3,4>)=\{\+l}A{i,,p,4>) : (20) 

the 'angular' part of the 'angular to 'radial' split, while the 'radial' coordinate x obeys, for our conformally-ordered free 
problem, 

-^—'{*t xmx <m\ + 5™*> + + Mn-l}{n-2 } _ 2^ 

sm d xcosx<ix I dx J snrx cos 2 x 2n - 3 h 2 

Now, (20) - analogous to the SHM part of the spherical harmonics equation - is a somewhat more complicated but still 
standard equation solved by the Wigner-D functions (see Appendix C). The quantum number I is interpreted as 'total 
isospin', i.e. in H-coordinates (Sec 1.2), the total angular momentum of the posts, and in K-coordinates (Sec 1.2) the 
picked-out 'axe blade' triangle subsystem. The I3 and Y quantum numbers are most candidly as the eigenvalues in 

X3DOA, fi, 0) - I 3 D(V, P, 0) and 5>D(V, 0, 4>) = X D (^, (3, 0) . (22) 

I 3 is thus interpreted as '3-componcnt of isospin', i.e. in H-coordinates the counter-rotation of the two posts relative to 
the crossbar, and in K-coordinates, the co-rotation of the face of the 'axe blade' and the handle relative to the depth 
of the blade (see Fig 1.8 for nomenclature). Also, Y is thus interpreted as the 'hypercharge' i.e. 'extra angular charge' 
quantum number, i.e. in H-coordinates the co-rotation of the two posts relative to the crossbar, and in K-coordinates, 
the counter-rotation of the face of the 'axe blade' and the handle relative to the depth of the blade. The ranges for the 
quantum numbers are 21 £ No and 2l3,Y £ Z with —I < I3, Y/2 < I. Finally, the 'radial' x -ec l ua tiori - analogous to the 
associated Legendre equation part of the spherical harmonics equation - can be shown [37] to map to the hypergeometric 
equation. Thus it is solved by Jacobi polynomials (within which the Gegenbauer and Legendre polynomials are nested, 
see Appendix B). Wigner D-functions are a combination of elementary trig functions and, once again, Jacobi polynomials 
(see Appendix C); they arise also in the elementary study of finite rotations in QM. 



5 Free problem's general solution and quadrilateralland interpretation 

For useful contrast, specializing (86) and (88) to the N = Y — 2 = 1 of pure-shape triangleland, 

£(k,l) = Ak{k + 1} , k £ N . (23) 

This is very familiar as a proportionality when k is denoted by 1 or J in the rigid rotor (though for triangleland itself the 
quantum number is denoted by S for 'shape'). The degeneracies are 

V(k,l) = 2k + 1 (24) 

which is also very familiar [SU(2) multiplets]. Thus this model has 1 ground state 's-orbital', 3 first excited state l p- 
orbitals' and 5 second excited state 'd-orbitals' using the spectoscopic notation familiar from Atomic Physics. For the 
triangle these are mathematically the same as for the atom but have the distinct physical interpretation provided in 
[15, 25]. Free 4-stop metroland has the above eigenvalues and multiplicities too [14, 25], only now the quantum number 
involved is denoted by D for 'dilational'. Free Y-stop metroland [14, 17] exhibits the reasonably well known 'rotor in Y— 1 
dimensions' pattern, which is also the mathematical basis of the simplest 'Y — 1-dimensional analogue of the periodic 
table' with Y - 1 'p-orbitals' and Y{Y - l}/2 - 1 'd-orbitals'. 
Next, specializing (86) and (88) to the N = Y - 2 = 2 of CP 2 [67, 37]), 

£{k,2) = 4k{k + 2} (25) 
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with degeneracies 

V(k,2) = {fc + l} 3 . (26) 

This has 1 ground state 's-orbital', 8 first excited state 'p-orbitals' and 27 second excited state 'rf-orbitals'. In this case, 
k := I + Y/2 + n, indeed motivated by being the sole functional dependence on the quantum numbers in the expression 
for the energy, with n <G No the degree in cos2x of the Jacobi polynomial in the corresponding expressions for the 
solutions at the energy level in question. This has some parallels with the principal quantum number of the atom and 
with its counterpart for the isotropic HO; these are all 'radial node counting' quantum numbers. The role of 'radius' in 
the current pure-shape problem is played by the \ coordinate. [I.e it is radial in contradistinction to the Euler angles' 
SU (2)-angularness, all to be taken within the context of the geometrically of configuration space rather than of the physics 
in space.] 

However, the above eigenvalues for Berger's and MacFarlane's treatments of CP 2 only involve Laplacian operator ordering 
and no factor of H 2 /2. Moreover, Ric(Q; M] is but constant here, so for quadrilateralland one can just take on MacFarlane's 
equation for a shifted energy as per above: 

E Uni = h 2 £/2 - h 2 2n{n - \}{n - 2}/{2n - 3} . (27) 

The second term here is —Ah 2 for the quadrilateral. Thus 

E Uni = 2h 2 {k + l} 2 , k e N . (28) 

Next, capitalizing on the split and identification at the end of Sec 4, the general wavefunctions for quadrilateralland are 

If {2{n + I + 1} + |Y|}r(n + 21 + |Y| + 2)n! T(I + Y/2 + l)r(I - Y/2 + 1) } 1/2 



tF i T(n + 21 + 2)r(n + |Y| + 1) T(I + I 3 + l)r(I - I 3 + 1) 

sin 2I X cosm xP (2i+i.|Y|) (cos 2 x )exp(i Y V/2)d Y/2 , l3 (/?)exp(i I 3 0) . (29) 

Key 22 This is obtained by interpreting MacFarlane's study [37] of the QM of CP 2 in quadrilateralland terms and also 
explicitly providing normalization coefficients for our convention (that he omitted) since we need them for our subsequent 
calculations. The mod bars come from the need for Re a, Re/3 > — 1 in the theory of the Jacobi polynomials (Appendix 
B); in this way we offer a minor correction of eq. (85) of MacFarlane [37]. The rest of the Y's present do not need mod 
bars by symmetry and the range of definition of the Wigncr-d and exponential functions; this is likewise for all of the Vs. 



For contrast, for the pure-shape versions of 4-stop metroland and triangleland are, respectively, 

*Dd«y Dd (e,$) cxP^(cos0)exp(id^) , (30) 

* sj cxr S j(e,$)cxP!j | (cose)cxp( i j$) , (31) 

for P now the associated Legendre functions. Thus we have a pure-ratio factor Rat(0) and a pure- relative-angle factor 
Ang($) in the triangleland case. 3 This feature is indeed repeated in quadrilateralland: a pure-ratio factor Rat(x, /?) and a 
pure-relative-angle factor Ang(0, ip). Note that while the Wigner D-function is a useful solving package, it is not aligned 
with this useful interpretational split, though, since it is itself separable into functions of all three variables, there is no 
problem in refactorizing the wavefunction. 

Thus both 4-stop metroland and triangleland free pure-shape problems have rigid-rotor mathematics, but the 4-stop 
metroland case is physically a relative dilatator and the triangleland case is a mixed relative rotor relative dilatator. 
This is probably best called a 'rationator'. Quadrilateralland is again a rationator, and one no longer constrained to 
obey SU (2) mathematics so that it is not the same as a rigid rotor. The symmetrical, alias Lagrange, spinning top with 
1\ = h 7^ I3, does itself involve at the quantum level eigenfunctions based on the Jacobi Polynomials [68], but these are 
mathematically distinct from the current paper's in greater detail. 

Parallelling further 4-stop metroland and triangleland calculations in [14, 15], e.g. 

sin 2I X cosl Y l X Pi 2I+1 ' |Y|) (cos2 X )cos(Y^/2)d Y/2j3 (/3)cos(l30) (32) 

can be recast as 

f T Y /2(n 2 • n 3 )T Y/2 (ni • n 3 ) 1 f Ti 3 (n 2 • n 3 )Ti 3 (n x • n 3 ) 1 

{l-n 2 }V 3 'Pf +1 ' |Y|) (2n 2 -l)| -T Y/2 (n 2 • n 3 )r Y/2 ( ni • n 3 ) |d Y/2il3 (arcta n (n 2 / ni ))| +Ti 3 (n 2 • n 3 )7i 3 (m • n 3 ) J 

(33) 

for T P (X) := yjl - T p {X) 2 and T p (X) the Tchebychev polynomial of the first kind of de gree p in X (see Appendix B). 

3 Pure-ratio here means a purely non-angular ratio, angles themselves of course being expressible in terms of certain other ratios; this all 
refers to the status of these coordinates interpretations in space; e.g. 8 is an angle in configuration space, but is physically a function of a 
purely non-angular ratio in space itself. 
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6 Visualization and discussion of first few wavefunctions 

This account parallels [14] and [13, 15] of 4-stop metroland and triangleland respectively. There, immediately visualizable 
2-d tessellations were available. However, now for quadrilateralland, we are at some disadvantage. Though at least some 
of the simpler wavefunctions can be viewed without loss with some dimensions suppressed. 

The ground state is 

s = *oooo = Vl/ir , const. (34) 
Thus it favours no particular regions or directions, 4 here meaning types of quadrilateral. 
1) The first harmonics are the octet 

PA :=*iooo = {2A}{l-3cos 2 X } , (35) 
Ps := *oioo = {2/7r}sin 2 xcos^ , (36) 
:= *oi-io = {V2/n}sin 2 xsm /3 cxp(-i<j)) , p s+ := * 1 1 o = -{V2/7r}sin 2 x sin/3 exp(i (j>) , (37) 
Pn '■= a - A i — {v / 3/7r}sin2xcos§ exp(— i (j>/2)exp(iip/2) , 
p p :=f i i] = {v / 3/7r}sin2xcos| exp(z</)/2)exp(ii/'/2) , 
Ps_ ■= a _i -i = {v / 3/7r}sin2xsin^ cxp(— i (/>/2)exp(— iip/2) , 

= {\^/7r}sin 2 X sinf cxp(i 0/2)exp(-i V/2) . (38) 




Figure 1: Gell-Mann's eightfold way multiplet familiar from Particle Physics, s is strangeness and q is charge, p and n arc the proton and 
the neutron. 

Note how quadrilateralland's orbitals, unlike those of the atom, its SO(d) generalization in dimension d, TV-stop metroland 
[14] and triangleland [13, 15], are no longer simply related to an obvious Cartesian space's axes. Thus they require a 
different kind of nomenclature for their labels. Wc still prefer sine and cosine combinations to ±'s, 0/-'s and p/n's, in 
parallel with the preferred representations of the atomic orbitals. See Fig 1 for the labelling nomenclature used in (34-38). 

Alternative specifically-quadrilateralland names for these orbitals are given in Figure 2 which sketches each's p.d.f. 
over our complex-projective chopping board generalization of Kendall's spherical blackboard in the case for which the 
Gibbons-Pope type coordinates are adapted about a Jacobi K-tree. 

By this stage we can see that the p's can be identified octet of shape variables (modulo normalization). I.e. we have 
the 'theorem' 

p r cx s r . (39) 

This is modulo taking sine and cosine combinations of the exponentials in the expressions for n and p, for E± and for 
H± and re-ordering the basis to (E c , E s , E , (n,p) c , (n,p) s , 5 C , 5 S , A). It is along the same lines ('naming polynomial' [14]) 
as how using Cartesians causes atomic orbitals to be the functions they are named after (subject to minor conventions 
such as that d z i is a contraction of d 3z 2_ 1 ). It is clear then that the deep-seated way of labelling orbitals for QM based 
on SU(n) is in terms of homogeneous polynomials rather than Cartesian axes (though both coincide for the well-known 
SO(3)-SU(2) case). Thus we have obtained a distinct, specifically quadrilateralland-based nomenclature for the orbitals 
in terms of the shape quantities. All in all, the 8 Gcll-Mann quadratic forms are natural successors of triangleland's 
3 Dragt quantities (that E.A. also termed Pauli quadratic forms): both sets are shape quantities, good for kinematical 
quantization and then a natural choice of labels for the quantum states for the corresponding free problem. 

With respect to the SU(2) privileged in this presentation, we have a singlet (A), a triplet (the E's) and two doublets 
(p/n and the S's). The particular symmetries of the states are pa is SU(2) x U(l) symmetric, ps is U(l) x U(l) 
symmetric, the ps ± are U(l) symmetric and the rest have no continuous symmetries. 

4 This is an ubiquitous feature - c.f. the ground state that is constant over the (k-)sphere for standard rotors and the corresponding .W-stop 
metroland and triangleland problems. 
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a) X=0 b) 




Figure 2: Using Fig I.8's complcx-projective chopping board geometrical representation and its associated nomenclature, a) pa's is 'long-and- 
short handled', though it slightly favours long over short by its nodal curve lying slightly into the short hemi-CP 2 . b) p-£ 's is 'pick-or-executioner 
bladcd'. c) ps ± and sinc/cosinc combinations thereof are 'regular'; further distinctions between these and within each of the below depend on 
additional relative-angle information, d) p^ Q or _ and sine/cosine combinations thereof are 'moderately pick-or-executioner' e) p n and pp and 
sine/cosinc combinations are 'all bar moderately pick-or-cxccutioncr'. x2 denotes that there are 2 distinct states with that pure-ratio profile. 
See Chapters 14 and 15 of [25] for 3- and 4-stop metroland and triangleland spherical blackboard equivalents of this Figure. 

7 Extension to scaled quadrilateralland QM (Key 23) 

Firstly, one gets a much closer match to GR quantum cosmology if one extends to the scaled C(CP 2 ) as in [1]. This is 
according to the following correspondences: configuration space radius p to scale factor a, shape degrees of freedom to 
inhomogencities, 'energy equation divided by moment of inertia' to 'Fricdmann equation post use of energy-momentum 
conservation equation' (see Chapter 5 of [25] for details). The RPM potential can furthermore be chosen so that this 
analogue Fricdmann equation parallels quantum cosmological scale dynamics. This has the further qualitative benefits 
that the associated small inhomogencity mathematics is a lot more tractable for these models than in the actual Halliwelh- 
Hawking scheme for GR itself. 5 

7.1 Kinematical quantization 

Lemma. Given a kinematical quantization algebra C for a shape space, then if the corresponding relational space has 
no 'extra' symmetries, the kinematical quantization algebra of the corresponding relational space is C © CXff for dff the 
'radial'/R+ problem's affine algebra. 

Example 1) Scaled iV-stop metroland is exceptional due to possessing a number of extra symmetries. Moreover, this 
particular case's mathematics is, of course well-known (c.f. Sec 1.2): the totality of the cone's symmetries that are not 
shape space symmetries are the translations. E.g. for 4-stop metroland, one has G can (S) V* = {50(3) (S) M 3 } (S) M 3 , which 
is mathematically a Heisenberg group. This has a second SO (3)-vector commutator and the standard commutation 
relation between the two conjugate vectors. 

Example 2) Scaled triangleland is also exceptional, working mathematically just like the n — 3 case of the preceding but 
physically the p l and iTi are now, rather Dra r and Iip ra . 

Example 3) Finally, for iV-a-gonland, N > 3, by Sec 2 and the Lemma, we take G can (S) V* = SU(n) (S) IHP(C", 2) (S) dff. 
Here, the sole extra nontrivial commutation relation is the affine one, [p, tt] = ihp. [n has to be represented as —ihpd p in 
order to succeed in being self-adjoint.] 

7.2 Conformal ordering and the TISE 

A^ (cp „_ 1} = A c(CP -i)f - 3{2n - 3}/V for A C(CP „- 1} = d 2 p + 2{n - Ijp' 1 + p- 2 A CP „-! . (40) 
Then the TISE is 

-{d 2 p + 2{n - l}/)- 1 + p- 2 {A CP „-! - 3{2n - 3}/4}* - 2{E Vni - V}/h 2 . (41) 

5 The classical equations for this arc in [25], at least for p, Z coordinates. The C(CP 2 ) equations in terms of p and Gibbons-Pope type 
coordinates may be new to us. We do not claim more than the radial conformal Killing vector, though it is a gap in our understanding as to 
whether this case furnishes more. 
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7.3 Scale-shape separation of the TISE 



Firstly note that the coning construction by which RPM's incorporate scale does not really care about the nature of 
the shape part, so the new split-out part comes out much the same as in [17]. For V = V(p) alone separability ensues 
= 5(shape alone) 7\L(p)]. 

p 2 K" + 2{n - 1} P W + {{2p 2 {£ Uni - V(p)}/h 2 - C}K = ; (42) 

as a new feature from quadrilateralland upwards, the separated-out shape part gives a different constant energy shifting, 
rather the conformal-ordered, pure-shape problem, 

{A c(CP „- 1)+ C-3{2n-3}/4}S = , (43) 

Comparing the first of these and the equation in Appendix A, we get that 

C = Ak{k + n- l} + 3{2n-3}/4 = 4fc{fc + 2} + 9/4 for quadrilateralland . (44) 

V = k constant gives but an equation that maps to (Appendix D) the Bessel equation, 

p 2 K" + 2{n - l}pK' + {2{E Vni - n}p 2 /h 2 - 4k{k + n - 1} - 3{2n - 3}/4}"K = , (45) 

whilst V = Ap 2 gives but an an equation that maps to (Appendix D) the associated Laguerre equation, 

p 2 K" + 2{n - l}pK' + {2{E Vni - K }p 2 /h 2 - 2Ap 4 - Ak{k + n - 1} - 3{2n - 3}/A}K = . (46) 

This is similar to the radial equation for each of the atom and the isotropic HO. 

7.4 General solution 

Thus the solution to the first of these for general A-a-gonland is 

K ex p{3„-2}/ 2j±v ______ ^2{E Vni - K }p/h) = K ex P {3 - 2}/2 J ±V 2{9+8fc+2fc2} (y/2{E Vai n}p/h ) (47) 

with the second equality specializing to quadrilateralland. This is simpler than the below, but has uncontained/un- 
normalizable character which is undesirable, so we mostly use the second example. 
For that, the general- N wavefunction is 

Rn(p) OC / ,{3-2n}/2+A„ L A„ {uJ p 2 / h)^^ p 2 /2) . (48) 

Then the quadrilateralland case is just the A„ to A subcase of this. Here, E — hiu{2{ii + A} + 1} for n e N, 
A„ := ^J{2n - 3}{2n - 9/4} + 16T{T + n - 1}} = ^45 + 64T{T + 2}/2 for the quadrilateral, and uj := VK = s/2A. 
These are roughly like the well-known radial profiles for atoms (themselves described by associated Laguerre functions) 
at least in terms of numbers of peaks and nodes. One qualitative difference is that atoms and A-stop metroland have the 
outer peak of the "2s" orbital p.d.f. much larger than the outer one. This corresponds to the most of the "2s" orbital 
lying outside the "Is" one, whereas this paper's A-a-gonland models have these two p.d.f. peaks of the same area as each 
other to within a few percent. 



8 A family of useful integrals (Key 24) 

We consider integrals of the form 

(Vi|0|V 2 ). (49) 

Subcases of these include the overlap integrals (for which the inserted operator O = id) and the expectation values (for 
which ipi = tp 2 )- Specific cases include (tpi | a n | ^2 ) (for powers of one's model's scale variable, a), and (tpi \ cos a | ^2 ) or 
(ipi I cos 2 a I ip2 ) for 01 a non-scale (i.e. preshape) variable. 
Three applications of these integrals are as follows. 

Application 1) expectation and spread of the scale and non-scale quantities in question. 

Application 2) Time-independent perturbation theory about e.g. free or HO-potential exact solutions. For quadrilat- 
eralland, however, it is not presently clear whether the extra indices that the Jacobi polynomials possess substantially 
complicate these calculations relative to those for the atom/triangleland by providing additional types of transition chan- 
nels/selection rules. Thus we do not yet know how to proceed to the particularly significant second-order case of this 
application, since this involves unequal quantum numbers on the 2 input wavefunctions as per the well-known general 
formula 6 

Ei il --E s , n , Es * Esi m v ' m 2 /{Ez ^> (50) 

6 Here, H and Q are multi-indices running over all of the system's quantum numbers. 



9 



for all that the simpler first-order formula 

E^ = (n\V'\Q) (51) 

is under control. We leave this point to a subsequent paper, noting that no such extra effects have for now been reported 
[53] for the intermediate-difficulty (Gcgenbauer polynomials, see Appendix B) problem of the Stark Effect in higher 
dimensions. 

Application 3) The really significant application for the current program is, moreover (by Sec 1.1) time-dependent per- 
turbation theory on the space of shapes with respect to the emergent time provided by the scale in the scale-shape split 
of scaled RPM models. This is useful due to its analogy with the Semiclassical Approach to the PoT and Quantum 
Cosmology. 

Note that both Applications 2) and 3) involve integrals of the specific form 

W|V>2> , (52) 

for V the perturbation part of the potential. 



Atomic Example 1) From the angular factors of the integrals trivially cancelling and orthogonality and recurrence relation 
properties of Lagucrre polynomials for the radial factors [66] , 

(nlmMnlm) = {3n 2 -l{l+l}}a /2 and A nlm r = ^{n 2 {n 2 + 2} - {1{1 + l}} 2 }a /2 , (53) 

where ao is the Bohr radius of the atom. One can then infer from this that 1) a minimal characteristic size is 3ao/2 for 
the ground state. 2) The radius and its spread both become large for large quantum numbers; moreover, for these, Brown 
showed that the classical orbits are well-approximated [70]. 

Atomic Example 2) (n' 1 m' | cos# sp | nlm) and (n' l' m' | cos 2 9 sp | n 1 m) are 3-Y integrals [64] (i.e. products of three 
spherical harmonics, Yjk, the radial parts of the integration now trivially cancelling). Here 'sp' denotes that the angles 
are taken in the the spatial sense that is common elsewhere than in this paper. Then the general case of 3-F integral 
is known, having been evaluated in terms of Wigner 3j symbols [64] . Integrals for the present Paper's specific cases of 
interest are furthermore provided case-by case in e.g. [69]. The first of these integrals occurs in the Stark effect [71] (±1 
selection rule) , and the second in the calculation underlying both Raman spectroscopy [72] and Pauling's analysis of the 
rotation of molecules within crystals [73] (±2 selection rule). For comparison with the below quadrilateralland working, 
the first of these integrals has as its nontrivial factor J^P™ (X)X~P™(X)dX and there is then a recurrence relation 

(97)by which XPj(X) can be turned into a linear combination of P^,,(X), whereupon orthonormality of the associated 
Legendre functions (96) can be applied to evaluate it. The second of these integrals then requires two uses of the same 
recurrence relation. 



RPM Example 1) For 4-stop metroland, the relevant shape integral or perturbed-potential integral is just the 1 — > D, m 
— > d, cos 2 # sp — > cos 2 # of Atomic Example 2) [14]. 

RPM Example 2) For triangleland, the relevant shape integral or perturbed-potential integral is just the 1 — » S, m — > j, 
cos6> sp — > cos 6 of Atomic Example 2) [15]. 

RPM Example 3) Parallels of Atomic Example 1) are given in [14, 15] and represent estimations of the RPM Bohr 
configuration space radius (or, for triangleland, Bohr moment of inertia) analogue of Atomic Physics' Bohr radius. One 
can furthermore view this as part of the Peaking Interpretation of Quantum Cosmology (see e.g. [74, 39, 25]), where the 
lack of universe-measurements rather constrains other means of 'interpreting QM'. 

RPM Example 4) New to the present paper, the relevant quadrilateralland shape integral selection rule for shape quantity 
s$ = cos 2% as the inserted operator is AI = AI3 = AY = 0, An = ±1 or 0. It more closely resembles the Stark Effect due 
to the ±1 part of its selection rule despite how the inserted term itself looks more like 4-stop metroland's (i.e. quadratic 
rather than linear, like for the Raman Effect). It is a case in which the model being an A-a-gon presides over the particle 
number N being 4, beacuse the Jacobi polynomials themselves are in cos 2%. Thus this serves as the basic-variable 
analogue of the Legendre variable, so a 'square' insertion is in fact a linear power in the basic 'Jacobi variable'. It differs 
by additionally allowing for the non-transition, 0. This feature is new to CP 2 , arising from the first RHS term in the 
recurrence relation (93) for the Jacobi polynomials, which is clearly zero for § p_1 /]R p /the p-d atom since a = j3 from the 
Gegenbauer polynomial specialization (99) downward. This example therefore unveils a number of good fortunes in the 
standard atomic version of these calculations, thus serving as a robustness test for the atom. The known surviving terms 
include the following one that is subsequently used in this paper, 

Y 2 - {21 + l} 2 

<* nIl3Y |cos2 X |* nIl3Y ) = {2n + 2I+|Y| + 1}{2n + 2I+|Y|+3} (54) 
and the first two nontrivial transition terms, 



/T 1 o it \ /n{n+|Y|}{n + 2I+l}{n + 2I+|Y| + l} 1 

(*„_! 113 y| cos 2 X |* nIl3 y) = y {2n + 2 I + |Y| + 2}{2n + 2I + |Y|} 2n + 2I+|Y| + l (55) 
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alongside the n — »■ n + 1 of this. These results and the subsequent perturbation theory applications can be viewed as 

further robustness tests for the mathematics and physics of the (arbitrary-dimensional) atom. 

See e.g. [71, 66] for atomic counterparts and [14, 15, 25] for 4-stop metroland and triangleland counterparts. 

RPM Example 5) The relevant quadrilateralland scale integrals are 

which can be interpreted as a Bohr configuration space radius, and 

(p 2 ) = {ft/w}{2n + 2A„ + 1} (for any n) (57) 

by the obvious factorization into scale and shape parts and recurrence relation (109). See Part III of File R for 4-stop 
metroland and triangleland counterparts. 

9 Application 1) Expectations and spreads 

One of us gave these for the 'Dragt' shape quantities for triangleland in [15] and, in collaboration with Franzen, for 4-stop 
metroland in [14]. For triangleland, these expectations came out to be zero. 

Comparison of mean angle (e.g. roughly from the expectation of cos 6) and mode angle (from graphs along the lines of 
those in [13]) reveals the mean to be larger than the mode, but by not quite as much as occurs radially in the atom. 
This reflects that this case's Gaussianity suppresses the mean-shifting tail more than the radial part of the atom's mere 
exponential does. Expectations and spreads of $ are just like for previous Sec as the 6-integrals trivially cancel in each 
case. 

For quadrilateralland, expectation of cos 2\ [proportional up to an additive constant to the simplest 'Gell-Mann quadratic 
form' shape quantity (1-71)] is given by eq (54) and its variance (i.e. measure of spread) by 

Var n ii 3Y (cos2x) = (* n ii 3 Y| cos 2 2x |* n ii 3 y) - (*nii 3 Y| cos2x |* n ii 3 y) 2 = 

n{n + 21 + |Y| + l}{n + 21 + l}{n + |Y|}{2n + 21 + |Y| + 3} 2 {2n + 21 + |Y| + 4} 
+{n + !}{n + 21 + |Y| + 2}{n + 21 + 2}{n + |Y| + !}{2n + 21 + |Y|}{2n + 21 + |Y| + l} 2 
X {2n + 21 + |Y|}{2n + 21 + |Y| + l} 2 {2n + 21 + |Y| + 2}{2n + 21 + |Y| + 3} 2 {2n + 21 + |Y| + 4} ( ] 

As limiting cases of particular interest, 1) the ground state has expectation - 1/3 and variance 2/9. 2) Expectation = 
— l/{2n + 3}{2n+ 1} — > l/4n 2 as n — > oo and variance — > 1/2 if I = = Y is kept throughout. 3) On the other hand, 
if n = I = Y is kept and n is sent to infinity, expectation goes to - 3/25 and variance to 196/625. 

Finally, Varp = 3.40fi/o; for the ground state of the scaled quadrilateral with isotropic HO potential - a particular case 
of confinedness controlled by the steepness of the well. 

10 Application 3) Emergent time-dependent perturbations in the Semi- 
classical Approach to Quantum Cosmology 

Following on from Sec 1.1, the specific r-presentation of Y-a-gonland unapproximated h and 1 equations are a Hamilton- 
Jacobi equation with quantum correction terms, 

{d h S} 2 -ihd h 2 S-2ihd h S( X \d h \ X ) -fi 2 {(x|a h 2 |x) + A : (7V,d)h- 1 ( x |a h |x>}-^h- 1 fc(Y,d)a h 5 

+h 2 h- 2 {c(N, d) - ( x |Ai|x» + 2V h (h) + 2( X \J(h, l a )|x) =2E Uni , (59) 
and what is for now a fluctuation equation 

{1 - P x }{ - 2thd h \ X )d h S - h 2 {d h 2 \ X ) + k(N, cOh-^lx) + h- 2 A ; }|x) + 2{Vi(l a ) + J(h, l a )}|x)} = . (60) 

These equations result from those in [75, 48] via various specializations in e.g. [40, 20, 25, 27]. 

Then the move sketched in eq 1.14 turns the first equation into a QM-corrected form of the classical energy equation, 

{*h} 2 - 2ih*h( X \d h \ X ) - h 2 {{ X \d h 2 \x) +k(N,d)h- 1 ( X \d h \x)} -*fih- 1 fc(7V,d)*h 

+h 2 h- 2 {c(N, d) - ( x |Ai|x» + 2V h (h) + 2(x|J(h, l a )|x) =2E Uni , (61) 
and the second equation into a QM-corrected TDSE, the core of which is 
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ftd|x)/dt om(WKB) = - - 2 ,_, w ^ mN Ai|x) + Ah 2 (t cm < WKB ))| x ) . (62) 



h 2 (("(WKB)) 

[Some omitted correction terms however cause this to depart from being a TDSE. 



This is a toy model of the GR Tomonaga-Schwinger equation (given here in relational formulation, i.e. in terms of frame 
and not shift 

ih{5/5t™ - {SF»/St cm }M^}\ X ) = i?i GR |x) • (63) 

(additionally coupled to the quantum momentum constraint equation M.^\x) — due to no prior explicit classical 
reduction being known in this case). Such an equation has been considered in more detail in the quantum-cosmological 
setting by Halliwell and Hawking [48]. The minisuperspace counterpart involves but partial derivatives, no correction 
term F^M^ and no coupled equation. 

The |x) separates into a new t-oi-p = h part 1Z and the same shape part as in the first half of this paper. Thus we identify 
|x) =R|nII 3 Y). Furthermore N.B. that (x|0|x) involves integration solely over the /-space = S(N,d). Thus we only need 
our pure-shape useful integrals for this application. 

We now take the classical t, solve (62), and then re- investigate the h-equation with this approximate knowledge of 
|x) so as to allow the 1-subsystem the opportunity to contribute to a corrected emergent timestandard. This is of 
course along the lines of the ephemeris time procedures outlined in Sec 1.14. STLRC is based on giving everything the 
opportunity to contribute but then ditching contributions that turn out to be negligible to the currently requisite accuracy. 
This means quantum emergent time has to be different from classical emergent time in principle, since the former has 
different/additional quantum changes contributing to it. [1, 46, 27, 25] lay this out. The part of it we consider in the 
present paper involves integrating up the h-equation to obtain 

f m(WKB) = f 2 6h/-B ± - 4C for (64) 

B = -ih{2( X \d h \x) +h- 1 k(N,d)} , C ^ -2{W h - { X \J\x)} + h^h- 1 k(N,d){ X \b\\x) - {x\dl\x) +^ 2 cXN,d)} . (65) 

So what were ± pairs of solutions to the Hamilton-Jacobi equation 1.14 at the classical level are turned into more distinct 
complex pairs, which splitting is mediated by operator-ordering and expectation contributions to first order in ft. One 
also sees that the second-order contributions are another expectation, another ordering term and one that has one factor's 
worth of each. 

Key 25 This is of the general Machian form [compare the classical counterpart (1.109)] 

f m(WKB) =^[h,l,dh,| X (h,l))]. (66) 
Expanding out and keeping up to 1 power of ft, 

,om(WKB) ,cm(WKB) 1 f (x|J|x) «ft f dh f k(N,d) \\,nf^\ (d7\ 

t =t (0) + ^=J -^ dh ~ 4/ W- h {~^ + {xldhlx) \ + 0( h (6?) 

i.e., with comparison with the classical counterpart (1. 110) an 'expectation of interaction' (J) term in place of an interaction 
term J, and an operator-ordering term and an expectation term in place of a classical 1-change term. 

The simplest case of interaction potential is J = Ap 2 cos2x- For this, the ordering term comes out as, in the ^4 = 
case — iH{{n — l}/2£uni}lnh = — {zft/i^uni}lnh for quadrilateralland. For A > 0, it comes out as 

— ih{{n— l}/2i?uni}hi(h/ \J E — Ah 2 = -{zft/i?uni}ln(h/ \J E — Ah 2 for quadrilateralland. N.B. this term does not involve 
any kind of coupling to the 1-equation, unlike the next three terms considered. 
Via (54), the J expectation term comes out as 

If y 2 -{2/+l} 2 I 1 (B) \ VAh , r-r-, ,,1 , . 

27li{2n + 2I + |Y| + l }}{ 2n + 2I + |Y| + 3}} 72 U}{^ Uni -Ah 2 ' ^^^^] ^ 

B/A = {K 2 — Ki}/{Ki + K 2 }, so this factor is an approximate contents homogeneity smallness. 
The other first-order expectation term, (dh), is 

ih f dli dTZ dt cm < WKB ) 1 f dh mttmrn ^, r , , „ l4l 

"2 J ^' dr(WKB) dh ='^=2 J ^K*n^k { k + n-l }+ A^ 

= -{^Hk + n-l}A + Z^ + En ^ + J^ arctan ( ^h \ 
AV2E Vni VE Vni -Ah 2 A 4V2AV 2 \^E Uni -Ah 2 ) 
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The higher-order derivative counterpart of the preceding, which occurs to second order is also analytically computiblc, 
coming out as proportional to 

1 -3A£3 ni h 6 + {9E* ni + AAA n {AAA n + 3£ 2 ni }}h 4 - 8A£ Uni A 2 h 2 - 2A 2 n E 2 ni 3E Vll 




, arctan 

6^2 i?f} ni hVSu„i - Ah 2 2V2 

f2A„ AA n ETj ni 2AA n f AA n \ 2 \ 

Note: that the shape scale TISE possesses exact solutions of what the semiclassical approach's TDSE merely approximates, 
so Sees 4-7 will eventually furnish tests for whether the main, and then smaller, regime choices done in the Semiclassical 
Approach arc consistent. 

Also note that there is a widespread prejudice in Semiclassical Quantum Cosmology that expectation terms are always 
small. However, counterexamples to this were given in [20] for trianglcland (expectation of 1-Laplacian term when the 
wavefunction is an eigenvector of the Laplacian) . In the present Paper we have shown furthermore that in some regions 
of configuration space the expectations of the interaction J, of <9h and of d 2 also blow up, so that the negligibility of 
expectation terms suffers from a Global Problem. Moreover, having to keep average terms spells the ends to any claims of 
Semiclassical Quantum Cosmology being an analytically-tractable subject. It would then not only be a numerical subject 
due to its many-term equations but also be of integro-differential form reminiscent of the Hartree-Fock [76] approximate 
formulation of Atomic and Molecular Physics. In this case, around 1930 or so, the inclusion of expectation terms was 
found to be highly necessary in order to at all accurately reproduce atomic and molecular spectra observations from one's 
Quantum Theory. See [47] for more. 

In terms of emergent 'rectified time' t cm ( rcc ) g i vcn by d/dt em ( rec ) : = h 2 d / dt om ( WKB ) OTj in shorthand, © := h 2 *, the 
1-equation is now cleaner, its t-dependence now being in line with basic Physics' TDSE: 

^| X )/9t cm(rcc) - -{fi 2 /2}A!| x ) + V(*° m(rcc) )lx> • (71) 

It then makes sense [27] to recast the h-equation in this same time (recollect the classical-level motivation for £ em ( JBB ) : 
the simplifier of equations of motion. It just now turns out that QM implies a conformally-related rectified time is to be 
used at the quantum level instead. The h-equation is then 

{®\nh} 2 - 2ih( X \*®\x) - ft 2 {(xl* 2 |x> + k(N, d)< X |*|x>} 
-ihk(N,d)®lnh + h 2 {k(£) - (Ai) = 2{£ rec - V£ c - ( X \Vr\x) - (x\J rcc \x)} ■ (72) 
for * := © - ®lc\ and * := */© lnh(t om ( roc )). Integrating this gives 

f m(rcc) = J 2dh/h 2 {-B ± VB 2 4C} . (73) 

Expanding out and keeping up to 1 power of ft, 

,om(roc) ,om(roc) 1 f (x\J\x) „ ^ f dh (k(N,d) \ 2 

t =*«» + W J ^dh- Tj / +2(x\d,\x)\+0(H). (74) 

The ordering term to the rectified time then comes out as, in the A — case, ih{n - l}/2.E Uni li 2 = iH/E Uni h 2 for 
quadrilateralland. For A > 0, it is ih{{n - l}/2£ Uni }{l/2h 2 - {A/ £ Uni }ln(h/ \J 'E - Ah 2 )} = {^/£ Uni }{l/2h 2 - 
{A/E\jni}ln(h/\/ E — Ah 2 )} for quadrilateralland. The (J) correction term to the rectified time is, via (54) again, 



Y 2 - {21 + l} 2 1 f B 



(75) 



2^2 l{2n + 2I+|Y| + l}}{2n + 2I+|Y|+3}i \E Vni j ^urn - Ah 2 ' 
The other first-order expectation term, (d^), is now 

E 2 ^ 2 + 2h 2 k { k + n-^ { 2Ahy E ^ + ^ arctan ( 7 J^L=) . (76) 



2V2E 2 Di hVE- Ah 2 2V2A \y/E Voi -Ah- 

Its higher-order derivative counterpart that occurs to second order is now proportional to 

const + A{5E* Di + 16AA 2 + 20AE 2 ni A»}h 6 - 2AE Vni A n {ZE 2 ni + 4AA n }h 4 2AE 2 ni A 2 n h 2 A 2 E^ 

' ~ ^unih'V^Uni-Ah 2 
\[A ( VAh \ \ A„ AA n ( AA n 1 A AA 2A . , 2N ) , x 
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More advanced cases of coupled h- and 1-equation schemes are considered in [20, 25, 27, 20], albeit just for 3-stop metroland 
and a few trianglcland workings so far. 

N.B. the (J), (<9 n ), (<9 2 ) terms are all contributions to/mechanisms for backreaction. Morevover, the (J) integral backre- 
action mechanism allowed for quadrilateralland is forbidden by symmetry /selection rules in the case of the triangle. This 
triangle to quadrilateral difference directly reflects the selection rule that is afforded by the Jacobi polynomials but not 
by their (Gegenbauer and) Legendre specializations. 



11 Naive Schrodinger Interpretation 

Armed with the present Paper's wavefunctions, let us now complete the Naive Schrodinger Interpretation [77, 78] evalu- 
ation of 



Prob (Region R) oc / |#„ii 3 y| dfi . (78) 



R 



as set up in Paper I and now labelled with this paper's quartet of quantum numbers. 
Example 1) 

Prob(e-collincar) oc / |* n i l3Y | 2 d^ (79) 



C e 



/■e,7T— e,27T — e,37r — e,4-7r /•€,7r— e,2-7T pir 

Prob(e-collincar) oc / III l*mi 3 y(X; Pi 4>i V0| 2sm3 X cos \ dx sin/3 d/3d4>dip 

Jib=0,Tr-e,2TT~e.3Tr-eATT-e J <h=0,Tr-e,2iT-e J 0=0 iv=0 



for the C e supplied in eq 1.145. Thus, taking the Gibbons-Pope type coordinates of this series of papers, 

/>e,7T— e,27T — €,3tt — e,47r /•€,7r— e,2-7T pir 

'■0=O,7T-e,2ir-e,37r-e,4Tr-e J <p=0,Tr-e,2iT-e J 0=0 J x- 

oc e 2 i.e. the cross-section of the '2-lune' for all at-least U(l) x £7(1) symmetric wavefunctions, i.e. Y = = I3 . (80) 

This concerns highly non-uniform states according to the demo(4) measure of uniformity (c.f. Sec 1.5.3). On the other 
hand, one of the first Y = 0, I3 7^ solutions gives the additional factor {1 — 2/|e} i.e. a small decrease. This shows 
that adding 'isospin' has the same-sign effect as adding relative dilational momentum to 4-stop metroland or mixed shape 
momentum to triangleland [25] . 
Example 2) Using the T e supplied in Appendix LA, 



P / 1 4-7T p27T P 7T P7l/2 

Prob(e-closc to a +43 triangle) OC / | 1 1 3 Y | 2 df2 = / / / / I *n 1 1 3 Y (X> A </>) V0| 2sin X cos X dx SU1 /3 d/3 d(j) dip 

JT t J4>=0 J<j>=0 J/3=Tr-e J x =0 



oc e i.e. the width of the belt for all at-least SU(2) x U(l) symmetric wavefunctions, i.e. Y = = I = I3. (81) 
Example 3) Using the S e supplied in Appendix LA, 



Prob(e-close to the 1243-labelled square) oc / l^niia y| 2 dri 



/ / / / |* n „ 3Y (x,/3,0,V')| 2 sin 3 xcosxdxsin/3d/3d0d ? A (82) 

J^ =T r_ e J0=O,27r-e J f3=TT/2-a J x=T /4-e 

oc e 4 i.e. the size of the '4-box' for the ground state n = = I = I3 = Y. (83) 

This one is a question of maximal uniformity, in fact more sharply defined than the demo(4) measure can provide, since 
the squares are not the only configurations that maximize that. Note that these are not substantially more probable than 
for other regions of configuration space of the same size, so there is not a big peak on high uniformity, unlike in one form 
of Barbour's conjecture [39]. This concurs with Sec 6's analysis. 

For comparison, GR Cosmology Naive Schrodinger Interpretation calculations can be found e.g. in [77]. 



12 Conclusion 

12.1 Quantization via use of geometrical methods 

We provided kincmatical quantization of the relational quadrilateral. The pure-shape version is based on the quadrilateral 
shape space's CP 2 's isometry group SU(3)/Z 3 and the linear space of Gell-Mann quadratic forms provided by Paper I. 
This is a coherent extension of trianglcland's kincmatical quantization, but only once one has taken into account the M. 3 
vector to Pauli matrix map by which trianglcland's shape quantities form Sec 1. 6's space IHP(C 2 , 2). Using C n for the 
linear space is geometrically natural too, but leads to absolutist rather than relational physics. For triangleland, 
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IHP(C 2 , 2) = R 3 has more minimal dimension than C 2 , whilst for quadrilatcralland, this minimality is reversed. (Isham 
[57] suggested, but did not oblige, dimensional minimality for the linear space to be involved in kinematical quantization.) 
The scaled version of kinematical quantization is based on the extension of this to the cone over CP 2 . 

We then provided time-independent Schrodinger equations for pure-shape and scaled quadrilatcralland. The first of 
these can build upon, after quantum-cosmological conformal-term adjustment, MacFarlane's work [37]. In the free case, 
it separates in Gibbons-Pope type coordinates. The second of these shape-scale splits into a different adjustment of the 
preceding pure-shape problem and a scale equation that just maps to the Bessel and associated Laguerre equations for 
the free and isotropic HO problems respectively. 

12.2 Interplay with Atomic/Molecular Physics, Particle Physics and Shape Geometry 

Next, we considered useful integrals - the analogue of those used in the atomic Stark Effect and for expectations and 
spreads of radial and shape operators in Atomic and Molecular Physics. In the quadrilateralland context, these can 
be used for time-independent perturbations about exact solutions, and time-dependent perturbations as useful in the 
Semiclassical Approach to the PoT and Quantum Cosmology. Furthermore, expectations and spreads of scale and shape 
operators are now a concrete means of carrying out the Peaking Interpretation of Quantum Cosmology [74, 39, 25]. 

The QM of CP 2 serves as a robustness test of the (fc-dimensional) rotor and atom problems as follows. 

A) this problem's orbitals can no longer be labelled by an obvious surrounding Cartesian space's axes (i.e. p x ,p yi p z , ... for 
the atom) ; instead, it has an octet of 'p-orbitals' that bear the same group-theoretic relations as in Gell-Mann's eightfold 
way in Particle Physics. Two things have happened here. 

1) The Cartesian space axes once again generalize in this way via equivalence to the Pauli matrices [SU(2) adjoint rep]. 

2) CP 2 has in excess of the degeneracy of orbitals possessed by that dimension's maximally symmetric shape space, S 
(e.g. by 8 to 5 for the p-orbitals). 

B) The ±1 selection rule problem now comes with a non-transition term (0 selection rule) absent from the /c-dimcnsional 
atom. This is a direct consequence of the Jacobi polynomials being more general than the Gegenbauer polynomials 
[compare (93) and (99)]. By it, the result concerning the second-order nature of the first perturbation terms for the rotor 
Stark effect does not carry over to CP 2 . 

We used the complex-projective chopping board of Paper I as a back-cloth for discussing the wavefunctions, in parallel 
to the use of spherical blackboards in [15, 14, 25], the triangleland case of which origninates in Kendall's work on Shape 
Statistics. What RPM's give back to that subject are questions concerning Geometrical Statistics and analogues of Shape 
Statistics itself for GR, QM and Quantum Gravity. 

Cones over complex projective spaces (and quotients of complex projective spaces and cones over those two) have 
featured in the String Theory literature e.g. as toy models of orbifolds [79]. Connections between these and Mechanics 
have long been pointed out by Atiyah. 

12.3 iV-a-gonland generalization of this Paper 

As the present Paper makes clear, quadrilateralland is far closer to the general TV-a-gon in terms of resultant mathematics, 
so the present paper is also the true gate to the general iV-a-gon. This SSec's considerations also lead to e.g. more general 
robustness studies along the lines of [56], large- N considerations such as Statistical Mechanics, a Shape Statistics approach 
[33, 36] to Records Theory [80] and study of the behaviour in the large-TV limit. 

Now the kinematical quantization involves Isom(CP fe ) = SU(k + l)/Zfc + i and IHP(C™, 2). As regards nonminimality, 
./V-a-gonland takes after quadrilateralland, since dim(Adj(5J7(n)) = n 2 — 1 > 2n = dim(C") for n > 2 (i.e. N > 3). 

The argument for conformal operator ordering is certainly general enough to hold for all ./V-a-gonlands. The resultant 
TISE is given in [25] in complex Fubini-Study coordinates. Some further work on this equation is given by MacFarlane 
in [81]. In particular, there continues to be an analogue of the radial \ part that separates out, and this continues to map 
to the hypergeometric equation and thus to give the Jacobi polynomials. The other separated-out part at this stage is 
a generalized Euler angle part consisting of the SU(N — 2) analogue of the 5/7(2) Wigner D-function. Here for now we 
do not know the extent to which the requisite analogy has been tabulated yet; this reason and making the new points 
without overly complicating the calculations and the presentation is why for now we stop at quadrilatcralland rather than 
solving for any JV-a-gon. MacFarlane also explicitly treats the whole TISE for CP 3 (i.e. for us, pentagonland) using the 
counterpart of the Gibbons-Pope type coordinates for this case. 

Also, via (86, 88) and Sec 2, 

(Energy of fcth eigenstate of iV-a-gonland) , E = 2h 2 {k{k + n - 2} + n{n - l}{n - 2}/{2n - 3}} , k e N , (84) 
with degeneracies 

D(k, N 2) = {N 2}{N - 2 + 2fc} | ^±^= ^ j . (85) 
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In particular, one has 1 ground state 's-orbital', {N - 2}N = n 2 — 1 first excited state 'p-orbitals' and {N - 1} 2 {N 2 - 
4}/4 second excited state 'd-orbitals'. Thus, interpreting MacFarlane's [81] in pseudo-atomic whole-universe terms, the 
6-d CP 3 pentagonland QM has 15 'p-orbitals' and 84 'd-orbitals', the 8-d CP 4 hexagonland QM has 24 'p-orbitals' and 
200 'd-orbitals', and that the 8-d CP 5 heptagonland QM has 35 ! p-orbitals' and 405 'd-orbitals'. 

As regards comparing HO potential terms and CP fe harmonics, CP fe 's first harmonic has dimension {fc + l} 2 — 1 i.e. 
the adjoint representation's i.e. that of the SU (fc + 1) itself. The HO's are pure-symmetric, so there are {k + 2}{fc + l}/2 
- 1 anisotropic modes among these, so they only cover part of the possible first harmonics. Contrast with the AT-stop 
metrolands for which these are all of the second harmonics [14]. The allness is due to lack of spatial antisymmetry in 
metroland - it is about what is omitted principally, which is the {k + l}fc/2 antisymmetric polynomials. This is because 
jth order harmonics are a basis for jth order polynomials but these can include antisymmetric polynomials for spatial 
dimension > 1 and these are not among the HO potentials. 

Our useful s$ = cos 2\ insertion %-integral has a clear A-a-gonland counterpart. The \ P ar t of this is no harder than 
in the present paper, though the remainder is less well-known for the analogues of the Wigncr D-functions. 

The scaled A-a-gonland's scaled part continues to separate out and give an equation of the same general form as for 
quadrilateralland (or, for that matter, any other RPM): eq (42). 

12.4 Problem of Time Applications in this Paper 

1) We considered quantum Kuchaf beables for the quadrilateral by aligning them with the kinematical quantization 
algebra, so that they are the 8 SU(3) generators and the 8 Gell-Mann quadratic forms. The most natural language for 
expressing all 16 of these at once is in terms of Gibbons-Pope type intrinsic coordinates. 

2) We considered the Machian version of Semiclassical Approach around the quantum Frozen Formalism Problem; in 
particular, we consider Machian correction terms to the zeroth approximation (itself not Machian) for the WKB time. 
Namely, we considered an operator-ordering term that can be treated decoupled from the quantum 1-physics of shapes 
and three types of backreaction terms that do require solving the quantum 1-TDSE. 

3) This paper's wavefunctions complement Paper I's characterization of physical propositions in terms of geometrically- 
simple regions of configuration space so as to be able to conclude Naive Schrodinger Interpretation calculations. This 
included consideration of quantum-cosmologically relvant questions concerning maximal and minimal uniformity. 
Paper III will contain each of Histories Theory and Records Theory for classical and quantum quadrilateralland. Paper IV 
will combine these with the present Paper's Machian Semiclassical Approach as a more advanced example of the program 
in [82, 83, 21, 25, 84]. This is a useful advance for the reasons given in the last paragraph of Paper I's Conclusion. It 
includes promoting the semiclassical quantum Kuchaf beables resolution to a quantum Dirac beables one. 

Acknowledgements: E.A.: I thank those close to me for being supportive of me whilst this work was done. Professors 
Don Page and Gary Gibbons for teaching me about CP 2 . Dr Julian Barbour for introducing me to RPMs. Mr Eduardo 
Serna for discussions. Professors Marc Lachieze-Rey, Malcolm MacCallum, Don Page, Reza Tavakol and Jeremy Butter- 
field for support with my career. This work was funded by a grant from the Foundational Questions Institute (FQXi) 
Fund, a donor-advised fund of the Silicon Valley Community Foundation on the basis of proposal FQXi-RFP3-1101 to 
the FQXi. Thanks also to Theiss Research and the CNRS for administering this grant. 

A The CP N eigenspectrum 

From Berger et al [85], 7 consideration of the eigenvalue problem for CP N (i.e. {A + £}u = so as to compare with our 
own convention for the free TISE), 

(fcth eigenvalue of CP N ) , £{k, N) = 4fc{N + k} , k e N , and (86) 

(Degeneracy of fcth cigenstate of CP N ) , V(k, N) ( N + k k f _ ( N t- f 1 )* ( 87 ) 

(from relating the cigenspaces to spaces of homogeneous polynomials whose dimension is elementarily computiblc), and 
which then simplifies to the more computationally useful form 

V(k, N) = N{N + 2fc} {(N + fc - 1)!/N! fc!} 2 . (88) 

We note that eq. 65 ii) of Macfarlane [81] has typos in it and should be replaced by 

dim(2, 0, 2) = X>(2, N) = N{N + 1} 2 {N + 4}/4 (89) 

(to convert between notations, our N is his n). 

7 This includes correcting a typo in the latter and making the following minor clarification. For k = 0, in the conceptual form, the reasoning 
is that ( N -i 1 ) is not a possible choosing process and therefore zero, and so D(0, N) = ( N o 1 ) — = 1 as indeed befits ground states. 
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B Jacobi polynomials 

The Jacobi polynomials [86, 87, 88, 89] P^\x) are the terminating series solutions that solve 

{1- x 2 }y" + {(3 - a- {a + (3 + 2}x}y' + n{n + a + (3 + l}y = , a, /3 > -1 , (90) 

which is the hypergeometric equation (under the map x = 1 — 2-q 2 ), this being the most general second-order linear o.d.e. in 
the complex plane to possess three simple poles. It includes the Gegenbauer alias ultraspherical polynomials as a special 
subcase (a = (3), with both the Legendre polynomials as a special subcase of that (a = = j3) and the Tchebychev 
polynomials of the first kind (a = — 1/2 = j3) as special subcases of that [86]. 
The Jacobi polynomials are standardized according to P^ a '^(l) = ( n n Q ). 
By recasting (90) in the Sturm-Liouville form 

{{1 - x} a+1 {l + xY +1 y'}' + n{n + a + (3 + 1}{1 - x} a {l + x} fj y = , (91) 

one can read off that the weight function is {1 — x} a {\ + x}^ . The orthonormality relation is then 

^(l-,H. + .)V(.^. 5 ^^«\. for R^B^-l. W 
We also need the following recurrence relation: 

2{n + l}{n + a + (3 + l}{2n + a + /3}P^f (x) = {2n + a + [3 + \}{{a 2 - p 2 } + (2n + a + /3){2n + a + (3 + 2}x}Y^\x) 

-2{n + a}{n + (3}{2n + a + (3 + 2}p£f } (x) . (93) 

For useful comparison: spherical harmonics and the (associated) Legendre equation 

The associated Legendre equation is 

{1 - X 2 }Y XX - 2XY X + {J{J + 1} - j 2 {l - X 2 }- 1 ^ = . (94) 

This comes from the 9 part of the spherical harmonics p.d.e. (X = cos#; the <p P ar t gives just SHM). 

It is solved by the associated Legendre functions Pj^(X) for J e No, j € Z, |j| < J. We use the standard convention 
that 

^) = {- 1 } j {l-^}^{^y^ 2 -l} J } . (95) 

by which 



is a complete set of orthonormal functions for X € [-1, 1]. We also require the recurrence relation [88, 86] 

XP a\ m = {J-ljl + n^W + iJ + ljll^lrW 

J V ' 2J + 1 ' K ' 

For useful comparison: ultraspherical harmonics and the Gegenbauer equation 

The Gegenbauer alias ultraspherical equation 

{1 - X 2 }Y XX - {2X + 1}XY X + J{J + 2\}Y = (98) 

is solved boundedly by the Gegenbauer Polynomials Cj(X; A). The {k > 3}-d ultraspherical harmonics equation arising 
as angular part of higher-(i problems is straightforwardly separable by the ansatz and change of variables into simple 
harmonic motion and a sequence of Gegenbauer problems. Normalization for these is provided in e.g. [86, 88]. The weight 
function is {1 — X 2 }^ 1 / 2 between equal-A Gegenbauer polynomials. These furthermore obey the recurrence relation 
[86, 88] 

XCAX; A) = {J + l}g J+ i(^;A) + {2A + J-l}C7 J - 1 (X;A) (gfl) 

2{J + Aj 

Tchebychev polynomials of the first kind 

The Tchebychev polynomials of the first kind T n (x) — cos(n arccos(x)) are the solutions of the Tchebychev equation 

{1 - x 2 }y xx - xy x + n 2 y = . (100) 
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C Wigner D-functions 

These are not usually separately tabulated or studied as special functions. This is because (see e.g. [68]) they can be 
expressed in terms of more basic special functions, which are tabulated and studied as special functions, according to the 
following relations [68]. 

The Wigner D-function D^ k (a, /?, 7) solves the equation 

{df 3 + cot/3 + shT 2 ^ 2 + d 2 - 2cos^<9 Q & r }Y~ + 1{1 + 1}Y = . (101) 

This separates into two SHM problems and a d-function of /3 alone that maps once again to the hypergeometric equation 
and thus gives Jacobi polynomials: 

D mkK ^ 7) = exp(im 7 )dU(^)exp(ika) , (102) 



dU(/3) := "?| I! = "! ! Bin- k f cos^f P<^ m+k W) . (103) 



{l + k}!{l-k}! 

The orthonormality relation for the Wigner D-functions is then 



^ Lo Lo Lo D - kl ( °' * 7)D - k2 K * l)Aa * ^ ^ = 2 V {U+m^liu-mJKl^mJI^-m;}!^- ^ ' 

(104) 



D Bessel functions and associated Laguerre polynomials 

These are two families of confluent hypergeometric functions. 
The Bessel equation of order p, 

v 2 w tVV + vw. v + {v 2 — p 2 }w = , (105) 
is solved by the Bessel functions. We denote Bessel functions of the first kind by J P (w). The family of equations 

2 - 1 {l-2a}xy. x + {a 2 + f3\k 2 x 2 ^ -p 2 }}y = (106) 



x y, 



XX 



map to the Bessel equation under the transformations w = x a y and v — kxP. The subcase of this with CX = 1/2, (3 = 1 
and p = 1 + 1/2 for 1 e N arc the well-known spherical Bessel functions [86] . 

The associated Laguerre polynomials are terminating functions for the confluent hypergeometric ode. What we need 
about them probably is as follows. The associated Laguerre equation 

xy, xx + {a + l-x}y, x +ny = (107) 

is solved by the associated Laguerre polynomials L"(z). They obey [86] the orthogonality relation 

/ x a c^p{-x)Vf j {x)Lf j ,{x)Ax = unless /3 = (108) 
Jo 

and the recurrence relation 

xL%{x) = {2/3 + a + l}L2(aO - {/3 + l}L a 0+1 (x) - {/? + ajhf^x) . (109) 
The 2-d quantum isotropic harmonic oscillator's radial equation for a particle of mass [i and oscillator frequency to, 

-{h 2 /2^}{R rr + R r /r + m 2 R/r 2 } + fiuj 2 r 2 R/2 = ER , (110) 
maps to the associated Laguerre equation under the asymptotically-motivated transformations 

R = {hx/iiLu} lnll/2 exp(-x/2)y(x) , x = [llot 2 jh . (Ill) 

Thus is solved by 

R cx rl m lexp( A1 o;r 2 /2fi)Ll m| (^r 2 /fi) (112) 
corresponding to the discrete energies E = {|m| + 2r + 1}Huj for radial quantum number r e No [90, 65]. 
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